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formulation and asymptotic analysis 
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Based on the asymptotic analysis of ordinary differential equations, we classify all spherically 
symmetric self-similar solutions to the Einstein equations which are asymptotically Friedmann at 
large distances and contain a perfect fluid with equation of state p = (7 — with < 7 < 
2/3. This corresponds to a "dark energy" fluid and the Friedmann solution is accelerated in this 
case due to anti-gravity. This extends the previous analysis of spherically symmetric self-similar 
solutions for fluids with positive pressure (7 > 1). However, in the latter case there is an additional 
parameter associated with the weak discontinuity at the sonic point and the solutions are only 
asymptotically "quasi- Friedmann" , in the sense that they exhibit an angle deficit at large distances. 
In the < 7 < 2/3 case, there is no sonic point and there exists a one-parameter family of 
solutions which are genuinely asymptotically Friedmann at large distances. We find eight classes 
of asymptotic behavior: Friedmann or quasi-Friedmann or quasi-static or constant-velocity at large 
distances, quasi-Friedmann or positive-mass singular or negative-mass singular at small distances, 
and quasi-Kantowski-Sachs at intermediate distances. The self-similar asymptotically quasi-static 
and quasi-Kantowski-Sachs solutions are analytically extendible and of great cosmological interest. 
We also investigate their conformal diagrams. The results of the present analysis are utilized in an 
accompanying paper to obtain and physically interpret numerical solutions. 

PACS numbers; 04.70.Bw, 95.36.+X, 97.60.Lf, 04.40.Nr, 04.25.Dm 



I. INTRODUCTION 

There is great interest in spherically symmetric self- 
similar solutions to Einstein's equations because of their 
numerous applications in astrophysics and cosmology [l| . 
Indeed there is now considerable evidence for the "simi- 
larity hypothesis", which postulates that spherically sym- 
metric solutions may naturally evolve to self-similar form 
in a variety of situations. The status of this hypothesis 
has been recently reviewed by Carr and Coley [3| . In view 
of this, it is clearly important to have as complete a clas- 
sification of spherically symmetric self-similar solutions 
as possible. 

Spherically symmetric self-similar solutions have the 
feature that all dimensionless quantities can be expressed 
in terms of z = r/t, where r and t are suitably chosen 
radial and time coordinates. If the source of the gravita- 
tional field is a single perfect fluid, Cahill and Taub Q 
have shown that the only barotropic equation of state 
compatible with the similarity assumption has the form 
p = (7 — 1)^ for some constant 7, where ^ and p are 
the energy density and the pressure, respectively. In this 
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case, for a given value of 7, the solutions are described 
by two parameters and, providing one restricts atten- 
tion to shock-free perfect fluids with positive pressure 
(1 < 7 < 2), they have been completely classified. This 
classification has been achieved using a combination of 
the "comoving" approach (in which the coordinates are 
adapted to the fluid 4-velocity) and the "homothetic" 
approach (in which the coordinates are adapted to the 
homothetic vector). These approaches have been used 
by Carr and Coley (3] (henceforth CC) and Goliath et 
al. , respectively, although a full understanding of the 
solutions requires that one combines them 

A key step in the CC analysis is the derivation of 
all possible asymptotic behaviors at large and small dis- 
tances For positive pressure, there are at least three 
kinds of behavior at large spatial distances (usually cor- 
responding to the limit z —^ 00): (1) asymptotically 
Friedmann (1-parameter); (2) asymptotically Kantowski- 
Sachs (1-parameter), though these are probably unphys- 
ical for 7 > 2/3); and (3) asymptotically quasi-static 
(2-parameter). There are also two families of solutions 
which exist only when 7 > 6/5: (4) asymptotically 
Minkowski at infinite z (1-parameter); and (5) asymptot- 
ically Minkowski at finite z but infinite physical distance 
(2-parameter). At small spatial distances, the solutions 
are of four kinds: they contain either (a) a black hole 
singularity or (b) a naked singularity at finite z (but zero 
physical distance) or they can be connected to z = via 
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a sonic point, in which case they are either (c) static or 
(d) represent a perturbation of the Friedmann solution. 

The complete family of 1 < 7 < 2 solutions can now 
be found by combining the five kinds of large-distance 
behavior and four kinds of small-distance behavior. The 
way in which one connects the large-distance and small- 
distance solutions depends crucially on whether or not 
the solution passes through a sonic point. If the solu- 
tions remain supersonic (or subsonic) everywhere, then 
the small- 2: behavior is uniquely determined by the large- 
z behavior. However, if there is a sonic point, the behav- 
ior of the solutions is much more complicated because 
the equations do not determine their behavior uniquely 
there, so there can be a discontinuity. Indeed, only a 
subset of solutions are "regular" at a sonic point in the 
sense that the pressure gradient is finite and they can 
be extended beyond it. Because of this feature, the fam- 
ily of solutions with a regular center and a regular sonic 
point have a band structure, with the solutions which arc 
analytic at the center and sonic point forming a discrete 
subset of these @, [l^l ■ This feature is very important 
for naked singularity formation and critical behavior in 
the gravitational collapse of a perfect fluid [ll|, • 

In this and the accompanying papers [Tsj (henceforth 
Paper II), we take a first step in extending the classifica- 
tion of CC to the negative pressure case (7 < 1), using 
a combination of numerical and analytical studies. In 
fact, part of the work is already done, because ref. 
does include the asymptotic analysis for 7 < 1. How- 
ever, there are some errors in that work (see Appendix 
A) and the full family of solutions has not yet been ana- 
lyzed. This paper will focus particularly on the case with 
< 7 < 2/3. Although this equation of state violates 
the strong energy condition and was little emphasized 
by CC, it may be very relevant for cosmolo gy - both in 
the early universe (when inflation occurred IJ] ) and at 
the current epoch [l5| (when acceleration may be driven 
by some form of "dark energy"). In both situations the 
matter model exhibits "anti-gravity" but the dominant, 
null and weak energy conditions still hold. 

As in the positive pressure case, it should be stressed 
that our classification does not cover imperfect or multi- 
ple fluids or solutions with shocks. The accelerated ex- 
pansion of the universe can also be realized by a scalar 
field with a flat potential (i.e. quintessence). The non- 
existence of self-similar black hole solutions embedded in 
an exact Friedmann universe for this case was already 
demonstrated in ref. and this is clearly complemen- 
tary to the present work. 

There are several key differences between the < 7 < 
2/3 and 2/3 < 7 < 2 cases. A purely formal difference 
is that the limiting values of z for large and small spa- 
tial distances in the flat Friedmann solution are reversed: 
large spatial distances now correspond to z — > and 
small ones to z — > 00. The other differences have more 
physical significance. The first is that there are no sound- 
waves, since the sound-speed Cg = ^ dp/d^ = 0-^/7 — 1 
is imaginary. This considerably simplifies the analysis. 



since there can be no discontinuities and solutions are 
analytic everywhere. The second is that there is no ex- 
act static solution, although there are still asymptoti- 
cally quasi-static solutions (a point which was missed in 
ref. Q). On the other hand, the Kantowski-Sachs and 
asymptotically quasi-Kantowski-Sachs solutions now be- 
come physical. Finally, the solutions which were asymp- 
totically Minkowski at large distances are now replaced 
with solutions which arc asymptotically singular at small 
distances. 

Despite these differences, many features of the CC 
classification still apply. In particular, there is still a 
1-parameter family of solutions asymptotic to the flat 
Friedmann model at large distances and these are of 
great physical interest since they may have cosmologi- 
cal applications. In the positive-pressure case, some of 
these solutions are supersonic everywhere and contain 
black holes which grow at the same rate as the particle 
horizon [itI . [18| . Others represent density perturbations 
in a Friedmann background which always maintain the 
same form relative to the particle horizon p^ . However, 
recently it was pointed out that none of these positive- 
pressure solutions are "properly" asymptotic Friedmann 
because they exhibit a solid angle deficit at infinity [20| . 
They may still be relevant to the real universe (since ob- 
servations may not preclude such an angle deficit) but 
it would be more accurate to describe them as "quasi- 
Friedmann" . 

In the < 7 < 2/3 case, we will show that there arc 
genuine asymptotically Friedmann solutions. We will an- 
alyze these solutions numerically and exploit these results 
in Paper II [l^ to interpret the solutions physically . The 
key feature is the existence of asymptotically Kantowski- 
Sachs and static solutions, both of which are extendible 
analytically. This leads to the possibility of cosmological 
black hole, wormhole and white hole solutions. 

The plan of this paper is as follows. In Section II, 
we present the basic field equations for spherically sym- 
metric self-similar spacetimes. We describe some exact 
solutions in Section HI and solutions which are asymp- 
totic to these In Section IV. We summarize our results 
and discuss their implications in Section V. 



II. BASIC EQUATIONS 

We consider a spherically symmetric spacctimc with 
the line element 

rfs2 ^ _g2*(t,r)^^2 ^ g2*(t,.)^^2 ^ ^^^^ ^Ydn\ (2.1) 

where dfJ^ = dO'^ + s\t? Odip'^ . We take the matter field 
to be a perfect fluid with energy-momentum tensor 



(2.2) 



where p and /i are the pressure and the energy density, 
respectively, is the 4- velocity of the fluid, and we use 
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units with c = 1. We will adopt comoving coordinates, 
so that the 4-velocity is 



dt' 



(2.3) 



The field equations can then be written in the following 
form: 





m 



-d 



'^,tR,r 
2* 



(2.4) 

(2.5) 

(2.6) 
(2.7) 
(2.8) 



e-^^i?,/), (2.9) 



where a comma denotes partial differentiation and G is 
the gravitational constant. The first two equations cor- 
respond to the energy-momentum conservation and the 
next two specify the Misner-Sharp mass to. Five of the 
above six equations are independent. Throughout this 
paper, we call the direction of increasing (decreasing) t 
the future (past) direction. 

A spacetime is self-similar if it admits a homothetic 
Killing vector which is defined by 

A5m- = S.^M'^, (2-10) 

where denotes the Lie derivative along ,f . Cahill 
and Taub [3] first investigated spherically symmetric self- 
similar solutions in which the homothetic Killing vector 
is neither parallel nor orthogonal to the fluid flow vector. 
They showed that - by a suitable coordinate transfor- 
mation - such solutions can be put into a form in which 
all dimensionless quantities are functions only of the di- 
mensionless variable z = r/t. In this case, ^'^ is given 
by 



d 
dxf^ 



d_ 

dr 



(2.11) 



The line element in a spherically symmetric self-similar 
spacetime can be written as 

ds^ = -e"^(-Ue + e2*(--)dr2 + r^S^{z)d9.\ (2.12) 

A hypersurface S of constant z is called a similarity sur- 
face and is generated by the homothetic Killing vector. 
The induced metric on S is 



dsi 



-(1 



V^)e^'^di^ + r^S'^dn^, (2.13) 

where V = |z|e*~* is the speed of the fluid relative to 
the surfaces of constant z. The area radius R = rS must 
be positive but S and r need not be. (Note that CC 
adopt a convention in which r is always positive.) A 
similarity surface is spacelikc for (y^ — l)e^* > 0, timelike 
for (1/2 - l)e2* < and nuU for (1 - V'^)e'^^ = 0. In the 



latter case, ^'^ is also null and the similarity surface is 
called a similarity horizon. 

It is helpful to rewrite Eq. (|2.12p in the conformally 
static form: 



-2* > (1 + V)dT + ^dz\ |(1 - V)dT - ^dz 



z^S^dQ^] 



(2.14) 



where the conformal factor depends only on t = In |t | . To 
obtain the causal structure of the self-similar spacetime, 
we can then use the analogy with the static spacetime. 
This is because the homothetic Killing vector, similarity 
surfaces and similarity horizons are the counterparts of 
the Killing vector, Killing surfaces and Killing horizons 
in the static case (see e.g. ref. [22|). 

The Einstein equations imply that p, fi and to must 
have the form 



S-kGp 
2Gm 



W{z) 
rM{z), 



(2.15) 

(2.16) 
(2.17) 



where we assume that the energy density is non-negative 
{W > 0). Then the Einstein and energy-momentum con- 
servation equations reduce to ordinary differential equa- 
tions for the non-dimensional functions with respect to 
the self-similar variable z: 



2P + P' 


= -(W + P)^', 




(2.18) 


-W 


= (W + P) (^*'4 




(2.19) 


M + M' 


^ WS^{S + S'), 




(2.20) 


-M' 


^ PS'^S', 




(2.21) 





= S" + S'-^'S' 


- *'(5 + 5') 


, (2.22) 


M 


= 5[1 + 2^6-2*5 


'-e-2*(5 4 


^ 5')(^,23) 



where a prime denotes a derivative with respect to In |z|. 

We assume the equation of state has the form p = 
(7 — 1)^, which is the only barotropic one compatible 
with the homothetic assumption. The dominant energy 
condition requires < 7 < 2. In this paper we exclude 
the value of 1 (corresponding to dust) since this needs 
special treatment and has been analyzed in ref. [2^ . We 
also exclude 7 = 0, corresponding to a cosmological 
constant, since this is incompatible with self-similarity. 
Equations (|2.18p and (|2.19p can then be integrated to 
give 

e* = coz2(''-i)/Tiy-(T-i)/T, (2.24) 
e* = ci5-2iy-i/^, (2.25) 

where cq and ci are integration constants. The velocity 
function V can be shown to be 



V 



Cl 
Co 



^(2-7)/7^-2yy(7-2)/7^ 



(2.26) 
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Equations (|2.18P " (|2.23p reduce to ordinary differential 
equations for S, S' and W: 



V^-{j-l)W' 7c?M^('^-2)/t 2(7-1) 



W 



2V- 



25*4 
S" 



S" = S' 



,(1-2 {^~l)W' 



1 



M = WS^{-iS' + S), 

M = s\l+ c-2^2(2-7)/7^2(7-l)/75'2 

-c^'^S^W'^'^{S + S'f 



(2.27) 

(2.28) 
(2.29) 

(2.30) 



From Eqs. (|2.29p and (|2.30p . we obtain the following re- 
lation between S, S' and W: 



WS'^{-lS' + S) = s 



1 + Co-2z2(2-^)/W2(^-i)/^5'' 



(2.31) 



Using this and Eq. (|2.29p . we obtain another constraint 
between M, S and W: 



M V 



7-1 + 



M \ 



1 - — I 
WS^J 



32) 



III. EXACT SOLUTIONS 
A. Friedmann solution 

The flat Friedmann solution in self-similar coordinates 
corresponds to 



(3.1) 



where 5*0 and Wq are constants. V is determined by 
Eq. ({2:26)) as 

Co 



and Eq. ((2?28)) holds trivially. From Eqs. (|2?27l) and 
(|2.3ip . 5*0 and Wq must satisfy 



,2u/('r-2)/7 



37 



U2 
^0 



25„4 ' 



(3.3) 



1 = c-25oX'^^(l-^) . (3.4) 



These equations are consistent provided 7 is not 2/3 or 
0. Since Sq and Wq are determined in terms of 7, cq and 
ci, where cq and ci are just gauge constants, there is no 
free parameter in this solution for a given value of 7. 

Introducing constants oq and 69 in place of Sq and Wq, 
we obtain: 



e = e = ao, 

I — (J 



M 



11-91 

4 



372a2 



463 



972a2|l-g|3' 



^2-39 



(3.5) 
(3.6) 

(3.7) 
(3.8) 
(3.9) 



Here the constants oq and can be chosen arbitrarily 
and 



37' 



(3.10) 



The constants cq and ci in Eqs. (|2.24p and (|2.25p are 
given in terms of oq and &o as 
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ao 



372^2 



(7-l)/7 



(1 — (7)2 ^3720, 



1/7 



(3.11) 
(3.12) 



where Eqs. (ISII1)-(IS!11), (EH) and are used. One can 

put the metric into a more familiar form, 



ds^ = -dP + P^dr^ + f^dn^), 
by introducing new coordinates 

— q\. 



(3.13) 



(3.14) 



It should be noted that r = 00 corresponds to r = 00 for 
the decelerating case (2/3 < 7 < 2) but to r = for the 
accelerating one (0 < 7 < 2/3). In the 7 ~ 2/3 case, the 
homothetic Killing vector is parallel to the fluid flow and 
these equations no longer apply [23 |. 

There exists a finite non-zero zi where V crosses 1. 
Figure [T] shows the conformal diagram of the flat Fried- 
mann solution for < 7 < 2/3, including the similarity 
surfaces. We can sec that the initial singularity is null. 
The similarity horizon z ~ zi corresponds to the cos- 
mological event horizon, while z — 00 and r = 00 is a 
regular center. The dotted and dot-dashed lines denote 
the singularity and infinity, respectively. 



B. Kantowski- Sachs solution 

The Kantowski-Sachs solution in self-similar coordi- 
nates corresponds to 



(3.15) 



5 



+0,t 




FIG. 1: The conformal diagram of the flat Friedmann solu- 
tion for < 7 < 2/3. There is a similarity horizon aX z = zi 
(> 0), corresponding to the cosmological event horizon, and 
a regular center at z = oo and r = oo. The initial singularity 
at t = is null, while null infinity is given by t = oo. The dot- 
ted and dot-dashed lines denote the singularity and infinity, 
respectively. The thin solid curves and lines denote similarity 
surfaces, i.e. orbits of z ^constant. 



Then Eq. (P?^ holds triviaUy and Eq. ([QT]) yields 



2-2/7 ^2 



OJ 



From Eq. (|2.26p . we have 
V = V^7?-^l\ 



Co 



(3.16) 



(3.17) 



and hence V ^ as z ^ oo. Eq. (|2.27[) is trivially 
satisfied to lowest order but determines Wq to the next 
lowest order: 



i-1 

7 



7^-1+2/7^ 



Equations (|3.16p and (|3.18p then give 
(2-37)(2-7), 



4(1-7) 



1, 



(3.18) 



(3.19) 



so we can obtain a physical solution only for < 7 < 2/3 
and there remains no free parameter in this case. This so- 
lution was first obtained in ref. [2I] , although that paper 
introduced different variables to deal with the unphysical 
solutions for the 2/3 < 7 < 2 case. 

The Kantowski-Sachs solution in more standard coor- 
dinates can be written as I21j 



2Gm 



= _(2-37)(2-7)^^.^^4(i^7)/7rfp 
7 

+ t^((i6'2 +sir^ 
^ 4(1-7) 

(2-37)(2-7)i2' 
4(l-7)2i 



(2-37)(2-7)' 



(3.20) 
(3.21) 

(3.22) 



V 



|^|^-(2-37)/7 

v/(2-37)(2-7)' 



(3.23) 



where f is a radial coordinate. As expected, V is only 
real for < 7 < 2/3. In the Kantowski-Sachs solution, 
the area of the 2-sphere with constant t and r docs not 
depend on r but only on i, so it expands with time. The 
topology of the constant t spacelike hypersurface is S"^ x 
R. 

There again exists a non-zero finite z\ where V crosses 
1. r = 00 or z = 00 for fixed t(> 0) can be analytically 
extended to negative r or negative z beyond z = ±00. 
Figure [2] shows the conformal diagram of the (extended) 
Kantowski-Sachs solution for < 7 < 2/3. The simi- 
larity surfaces are also shown. We can see that the ini- 
tial singularity is null. There are two similarity horizons 
z = ±zi, corresponding to two cosmological event hori- 
zons. The corresponding conformal diagram (Fig. 15) in 
ref. I22II is incorrect. 



-0,t 



-0,t 




FIG. 2: The conformal diagram of the Kantowski-Sachs solu- 
tion for < 7 < 2/3. The solution is analytically extendible 
beyond z — ±00 and r — ±co. There are two similarity 
horizons 2 — ±21, both corresponding to cosmological event 
horizons. The initial singularity at t = is null and future 
null infinity is given by t = 00. 



C. Absence of static solution 

The static solution would need to have 
S = So, Wo, 
in which case Eq. (|2.26p gives 



V = Voz^^-^y\ Vo = ^So^w^''-^^'\ 

Co 



Equations (|2.27p and (|2.3ip then reduce to 

•1-7 



= ^ ° + 2 



St 



7 



and 



(3.24) 



(3.25) 



(3.26) 
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respectively. For 7 > 1, Eq. (|3.25p can be satisfied and 
there is a static self-similar solution, but there is no such 
solution for < 7 < 1 However, as we will see later, the 
absence of an exact static solution does not preclude the 
possibility of an asymptotically static solution. 



IV. ASYMPTOTIC BEHAVIORS 

A. Friedmann asymptote for 2^0 

We now focus on self-similar solutions which are 
asymptotic to the flat Friedmann model at large spatial 
distances, i.e. in which <I>, 5", S, W and M approach 
the form given by Eqs. (|3.5p - ()3.9p as f — > 00. In this 
case, it is convenient to introduce new functions A and 
B, defined by 



S = SF{z)e^^'\ 



(4.1) 
(4.2) 



which describe the deviations from the flat Friedmann 
solution. Equations (P?^ and (P?^ can then be 

written as 



A' 



B" = 



7y^[37g^e-('^-^)^/'y - 2(g + 2B')] 
2[V^2-(7-l)] - 

-(B'-q) (l-q + B' + "^'^ 



A' 



V 



2B' + -A' ] (1 



!^^l'qg-2B + {l-2/y)A^ 
Oo 



7 

q + B') 



(4.3) 

(4.4) 
(4.5) 



respectively. Equation (|2.3ip yields another relation: 



97X11 
bo 



1 



1 



a'oil-qy 



,2-2g 



B' 



(4.6) 



The flat Friedmann solution is given hy A = B = for 
all z. Spacelikc infinity corresponds to z 00 for the de- 
celerating case (2/3 < 7 < 2) and z — > for the acceler- 
ating case (0 < 7 < 2/3). Note that V cx ^1-2/(37) _^ ^ 
at large distances if A and B are finite. This means that 
z = 00 (2/3 < 7 < 2) and z = (0 < 7 < 2/3) are 
horizontal lines in a conformal diagram. 

The asymptotic form of self-similar solutions which ap- 
proach the flat Friedmann at large distances can be found 
by neglecting the V~'^ term and linearizing the equations 
with respect to A, A' , B' and B" . (It is noted that B 



does not need to be small.) From Eqs. (|4.3p and (|4.4p . 
the linearized equations at spacelikc infinity become 



= 3g2(7-2)A + 4B' + -A', 



(4.7) 







B" + 'i{l-q)B' +[^—^]A' (4.8) 
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in both cases. These equations lead to two independent 
solutions and the general solution is given by a linear 
combination of these. 
The first solution is 

A{z) = Aiz2(2-37)/(37) + ...^ (4.9) 
B{z) = /3 + Biz2(2-37)/(37) + ...^ (4.10) 



where 



Bi 



(2-37)Si, 
ap^ 7(37 - 2) 
feo' 4(37 + 2) 



-2/3 



(4.11) 
(4.12) 



/? is a constant and the dots denote higher order terms. 
A and B converge at large distances in both the decel- 
erating (2/3 < 7 < 2) and accelerating (0 < 7 < 2/3) 
cases. This solution was first obtained in ref. [l3| for the 
radiation case (7 ~ 4/3) and in ref. [Toj for more gen- 
eral 7. However, it is not properly asymptotic to the flat 
Friedmann solution because Eq. (|4.2p shows that there is 
a residual e'^ term in S at infinity, so Eq. (|2.12p implies 
that there is a solid angle deficit We describe so- 

lutions with this asymptotic behavior as asymptotically 
quasi- Friedmann. 
The second solution is 



A 



A^z^-'--'^'^ 



1 



07 



(4.13) 
(4.14) 



where is a constant and x fa y means (x/y) — i- 1 in 
the relevant limit. The above solution was also found in 
ref. [2^. The decelerating solution must be discarded, 
since this diverges at large distances (z — s- 00). On 
the other hand, the accelerating solution is "properly" 
asymptotic to Friedmann since A and B both converge 
to zero at large distances (z — > 0), so there is no solid an- 
gle deficit. In the following, we generally omit the term 
"properly" and simply describe these solutions as asymp- 
totically Friedmann. 

We conclude that there is a 1-paramcter family of 
asymptotically Friedmann solutions at large distances for 
< 7 < 2/3. Note that, in the positive pressure case 
(7 > 1), it is well known that there is no physical solution 
which is exactly Friedmann at sufficiently large but finite 
distances [l3l- Although one might envisage attaching 
an interior black hole solution to an exterior Friedmann 
solution at a sonic point (since there can be a discontinu- 
ity there), no such solution is possible. This conclusion 
trivially applies when the pressure is negative (7 < 1) 
because there are no sound-waves in this case. 
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B. Friedmann asymptote for z ^ oo 

We are also interested in solutions which are asymp- 
totically Friedmann at small distances from the origin. 
We therefore seek solutions which have 



(4.15) 



for 2; ^ oo and < 7 < 2/3, where the constants 5*0 and 
Wo may be different from the exact flat Friedmann case. 
Equation (|2.26p implies 



Co 



(4.16) 



The exact Friedmann relation p.4p still applies but 
Eq. I|2.27p now yields to lowest order 



37 ° -f \W 



2S^ 



7c?Wo(^-')/'^ 
J4.17) 

Hence {W /W) — 2 is either proportional to 4/(37) 
or falls off even faster. Therefore Eq. (|2.28p implies the 
following asymptotic behavior: 



S 

w 



-2/(37) 



l + 5lz2(l-2/(37)) + . 



where 



Wnz' 



37 



1 + P^^^2(l-2/(37)) + . . . 



1-7 



1 



Wi 



2^4 



57(3 - 27) 



,(4.18) 
(4.19) 

(4.20) 
(4.21) 



All coefficients of higher order terms are determined from 
So and W^. Since the value for can be different from 
the exact Friedmann case, there is a 1-parameter family 
of asymptotically quasi- Friedmann solutions at small dis- 
tances. Since s-Oasz— i-oo,2 = ooisa vertical line in 
a conformal diagram. These solutions can be interpreted 
as self-similar models with a regular center as r — *■ for 
fixed t 7^ 0. On the other hand, in the limit t ^ for 
fixed r 7^ 0, they correspond to simultaneous big bang 
models in the comoving time-slicing. If Sq had the ex- 
act Friedmann value, we would have Friedmann itself, so 
this is the only solution which is properly asymptotically 
Friedmann. 



C. Static asymptote for 2 ^ 00 

As we have seen, there exists no static self-similar so- 
lution for < 7 < 2/3. However, we will show that there 
are still solutions which can be described as asymptoti- 
cally quasi-static. In this case, we assume 



WkWo>0, SKiSo>0 



as z ^ 00. Only if Sq and Wq are the same as in the static 
solution do we describe such solutions as asymptotically 
static. Equation (|2.26p implies 



V: 



y„^(2-7)/7, V, = £iz(2-7)/7Q-2w(7-2)/7^ 



Co 



(4.22) 

This solution can be compatible with Eq. (|2.27p if and 
only if 



-f W b 



(4.23) 



as 2 ^ 00. This combination appears in the second term 
on the right-hand side of Eq. (I2.28P and this can be re- 
garded higher order term. It then follows that 



S = So + Siz'-^/^ + 52^2(1-2/7) + . . . . 

From Eq. (|4.23p . W must be of the form 

W = Wo + W^izl-2/7 + 1^2^2(1-2/7) _^ . . 

where 

1 Wi 

7 Wo 

From Eq. (|2.3ip . we obtain the relation 

2 



4 

Do 



0. 



(4.24) 



(4.25) 



(4.26) 



WoSi = So 



1 



1 



■Wr 



2(1-2/7) C.2 







To lowest order Eq. (|2.27p becomes 



2V' 1 



IW2 
7 Wo 



.So 



i4Wo 



Si 



1-2/7 



c^iS'oWo 
(4.27) 



2/7 



2St 



and to the next lowest order Eq. (|2.28p becomes 



2(1-7) 
7 

(4.28) 



2(1--) 52 = - I 1 - - 

7 



1 - - j SiWi 

W2 



7W0 



i4.29) 



5i and Wi are determined from Eqs. and (|i?771) 

in terms of 5o and Wo 7 while 52 and W2 are determined 
by Eqs. and ^^M) in terms of 5o, 5i, Wo and Wi. 

All higher order terms are determined in terms of So and 
Wo- Hence there is a 2-parameter family of asymptoti- 
cally quasi-static solutions. Since there is no exact static 
solution, there are no asymptotically static ones. Also 

— > cx) as z — > 00 for < 7 < 2/3, so z = 00 is a 
horizontal line in a conformal diagram. 

If we consider the analytic continuation beyond z ~ 00 
into the negative z region, the metric should remain ana- 
lytic in terms of the local inertial Cartesian coordinates. 
For asymptotically quasi-static solutions, one should con- 
sider the analytic continuation beyond t = -1-0 because 
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S = R/r and W = SirGfir'^ are finite. For this purpose, 
it is convenient to see how the proper time t behaves for 
fixed r > 0. For this class of solutions, Eq. (|2.24p implies 



Coz-2(l-7)/7^a-7)/7^ 



(4.30) 



so T (X t 1+2/7. When we continue the solution analyti- 
cally to the negative t region beyond t = ±0, this means 
that T « iC^ltl-i+^/T, where a is a positive constant 
depending on r and the upper (lower) sign corresponds to 
the positive (negative) t. To express W and S as analytic 
functions of t, we use the following unique continuation: 

S = 5o±5i|2|i-2/^ + 52|z|2(i-2/7) + ...^ (4.31) 
W = Wo± W^i|z|i-2/7 + vi/2|z|2(i-2/7) + ...(.4.32) 

This expression gives a Taylor series expansion around 
r = in terms of r. Because of the presence of the odd 
powers of r, this is not reflection-symmetric about t — 0. 

On the other hand, in the limit r ^ oo for fixed t 0, 
one has a vacuum of infinite radius which is not asymp- 
totically flat because m/R approaches a non-zero con- 
stant. We call this quasi-static spacelike infinity. In 
Paper II, we will see that there are a class of solutions 
describing a Friedmann universe emergent from a white 
hole, where a Friedmann spacelikc infinity and a quasi- 
static spacelike infinity are connected. 



D. Kantowski-Sachs asymptote for z oo 

For the asymptotically quasi- Kantowski-Sachs solu- 
tions, we assume 



as z CO. From Eq. (|2.26p . wc have 



(4.33) 



V « V^z''-^'^ = ^-^S^^W^-^'^ z^'^'\ (4.34) 
Co 



and hence F ^ as z ^ oo. Equation (|2.27p is trivially 
satisfied to lowest order but to the next order requires 
either 



W 



2 CX ^2(3-2/7) 



(4.35) 



or that the right-hand side falls off even faster than this. 
Equations (P?^ and imply 

S = 5o^-l(l-f ^1^3-2/7 + 52^2(3-2/7) + ... 1^36) 

W = H/o^^(l + W^2z^'^"^/''^ + •••) (4.37) 
Equation ((2?3T|) yields 



(l-7)Wo53 = 5*0 



2w2(7-l)/7 g2 







1 + c-^'W, 



(4.38) 



while Eq. ((2?27| gives 



7 7X7/ 



2^4 



At second order Eq. (|2.28|) becomes 



(4.39) 



S2 



1 



{W2 + Si] 



(4.40) 

Hence the coefficients of all higher order terms are de- 
termined in terms of 5o and Wq. This means there is a 
2-parameter family of asymptotically quasi-Kantowski- 
Sachs solutions. If the values for S'o and Wq are the 
same as those for the exact Kantowski-Sachs solution, 
all coefficients of higher order terms vanish, so the solu- 
tion is also exactly Kantowski-Sachs. In other words, the 
only solution which is asymptotically Kantowski-Sachs is 
Kantowski-Sachs itself. On the other hand, even if the 
first order terms vanish, i.e. Si = 0, we can see that W2 
may not vanish, so S'o and Wq arc different from their 
Kantowski-Sachs values. Since 1/ ^ as z ^ 00 for 
0<7<2/3,z = (X)isa vertical line in a conformal 
diagram. 

For asymptotically quasi-Kantowski-Sachs solutions, 
one can consider the analytic continuation beyond r = 00 
because zS ~ R/t and W/z'^ = SnG^t^ are non-zero and 
finite for fixed t ^ Q. It is useful to see how the proper 
length A changes around r = 00 for fixed t ^ 0. Equa- 
tion (|2.25p implies 



CiZ 



(4.41) 



so A CX r3-2/7. This means that, when we continue 
the solution analytically to the negative r region beyond 
r = 00, one has A « ^Ct\r\^^~^^'^ , where Ct is a posi- 
tive constant depending on t and the upper (lower) sign 
corresponds to positive (negative) r. To obtain analytic 
expressions for W and S in terms of A, we use the fol- 
lowing continuation: 



S = Soz-\l±Si\z\^-^^^ + S: 



2\Z 



2(3-2/7) 



W 



Woz'(l + W^2Nl'*^"'/^^ + •••)• 



+ • -(4,42) 
(4.43) 



The above expression gives a Taylor series expansion 
around A 0. Since Wi = in the above expression, 
W/ z'^ = SirGt^ generically has an cxtrcmum at z = ±00, 
while zS = R/t does not. So this expression is not 
reflection-symmetric with respect to r ~ ±00. On the 
other hand, in the limit t —^ for fixed r 7^ 0, we have 
R and jj, 00, which corresponds to an initial 
singularity. In Paper II, we will see that this analytical 
continuation is crucial for obtaining black hole solutions 
embedded in a Friedmann background as well as a class of 
wormhole solutions connecting one Friedmann spacelike 
infinity to another quasi-Friedmann spacelikc infinity. 
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E. Constant-velocity asymptote for z ^ oo 

Here we seek solutions in which V tends to a finite 
positive value Voo as z — > cx). Differentiating Eq. (|2.26p 
and noting that V' /V — > as z oo, we have 



W 



1 



27 
i-7 



S' 



(4.44) 



Substituting this relation into Eq. (|2.28p . we find to low- 
est order 



5^ _ _f _ 4(1-7) 
5 ~ " 



7 



(2-7) 



2(2 -7^) ^ 
7(2-7) S- 



(4.45) 



To lowest order this ordinary differential equation has a 
general solution of the form 



where So^± are arbitrary constants and 



(4.46) 



P± 



-(2 - 7^) ± ^(l-7)(4- 87 + 472 -73) 



(6 - 57)7 

(4.47) 

One can show that the square root is real for < 7 < 
2/3. Also Eq. (|2.3ip is consistent only if the first term 
dominates the second term in Eq. (|4.46p . so Sq = is 
excluded. 

One can easily show that the first term on the right- 
hand side of Eq. (|2.27p converges to zero for < 7 < 2/3 
and so 



_ 7(1 - 7) + V(l - 7)(4 - 87 + 47^ - 7^) 
Voo- ^—^ . (4.48) 

This is plotted in Fig. [31 Note that the solution with 
similar asymptotic behavior for 7 > 6/5 [1] has the other 
branch of the square root. For this solution, we have 



(4.49) 



The condition V 
Wo: 



Voo yields a relation between Sq and 



Co 



(4.50) 



There is no other relation between them, so there is a 
1-parameter family of solutions with this asymptotic be- 
havior. Since V Voo < 1 as z — > 00 for < 7 < 2/3, 
z = 00 is a vertical line in a conformal diagram. 

We now discuss the physical significance of this asymp- 
totic solution. As 7 increases from to 2/3, p+ decreases 
from —1/2 to —3/4, while Voo decreases from 1 to 1/3. 
Because M oc z1+(6-57)p+/(2^7)^ we have R = rS ^ 00, 
fj. = W/{8T:Gr^) ^ 0, m = rM/{2G) 00 and 
2m /R = AZ/S* ^ 00 as ?' ^ 00 for fixed i 7^ 0. There- 
fore, this solution approaches a vacuum region with in- 
finite physical radius. However, because the fall-off of 



■ 




FIG. 3: Voo for the asymptotically constant- velocity solution 
as a function of 7. Although only the part of the solid curve 
with < 7 < 2/3 is relevant to the present analysis, it is 
extended to < 7 < 2 for completeness. The other root 
for Voo is also plotted with a dashed line but this may be 
unphysical. 



the energy density ^ is very slow, the solution is far 
from asymptotically flat. In fact, both m and m/R di- 
verge to infinity. This solution cannot be analytically 
extended beyond z = 00 because R diverges to infinity 
as r — > 00 for fixed t 7^ 0, while diverges to infinity as 
t ^ for fixed r 7^ 0. We describe these as asymptoti- 
cally constant-velocity solutions. The infinity reached as 
r — > 00 for fixed t ^ Q will be called the constant-velocity 
spacelike infinity (cf. the asymptotically Minkowski so- 
lutions for 7 > 6/5 described in CC). 



F. Singular asymptote for z 

We assume that InlS"! diverges as z z*, while V 
tends to a finite value • Then Eq. (|2.26p implies 



W 
W 



27 S' 



2-75 

and Eq. (|2.28p yields to lowest order 



S 

w 



5p|^|(2-7)/(6^57)^ 
^^|^|-27/(6-57)^ 



where 



Z = In- 



(4.51) 



(4.52) 
(4.53) 



(4.54) 



Equation f^^Tl} implies V;^ = 1, while Eq. (P?^ gives 

(4.55) 



6 — 57 
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where the upper (lower) sign corresponds to positive 
(negative) Z. From Eq. (|2.26p . we have 



%|2,p(2-7)/75-4^2(7-2)/7^^ 



(4.56) 



and hence jz*! is determined from Sq and Wo- To next 
lowest order, Eq. (|2.28p implies 

S = S'o|Z|(2-T)/(6-57) [l + 5-^|^|(2-37)/(6-57) ^ . . (\ 57) 



W = Wq\Z\ 

where 



-27/(6-57) 



l + M^l|Z|(2-3'')/(6-57) + ..(.] ^8) 



_ (2-7) 
Through Eq. (H^Sl), we obtain 

F=l + Fi|Z|(2-37)/(6-57) + . 

where 



and hence 



Vi 

Si = 
Wi = 



-25*1 



7 



(1-7) 
7(4 - 37) 



-Vi. 



4(2-7)(l-7) 
To next lowest order Eq. (|2.27|) then gives 

'6-57\ 7c2Wtr('-^)/^ 



V,=± 



■ 7 



2^4 



(4.59) 

(4.60) 
(4.61) 

(4.62) 
(4.63) 

(4.64) 



which is consistent with Eq. (|2.3ip . 

We describe these solutions as asymptotically singular. 
For fixed r 7^ 0, we have R = rS ^0, n = W/{87rGr^) -> 
00 and m = rM/{2G) w rMo/{2G) as Z ^ 0. From 
Eq. (|4.55p . we have 



m : 



rMa 
~2G 



S^^ ^oS'oR (4.65) 
6 — 07 



and so the mass is positive (negative) for Z ^ 0^ 
(Z 0~). (Note that S can be negative but R = rS 
is positive even after the extension to the negative r re- 
gion.) Since the physical properties of these two lim- 
its are very different, we distinguish these two cases as 
asymptotically positive-mass singular and asymptotically 
negative-mass singular. For the positive-mass singular 
case, Z ^ 0+, y ^ 1+ and m > 0; for the negative case, 
Z ^ 0-,V ^ 1- and m < 0. Although F ^ 1 in the 
limit Z ^ 0, the first component of the metric (|2.13p is 



2* 



2cl\z, 



|4(7-l)/7^2(l-7)/7 



-(2-7)7(6-57/ 



f4.66) 



which tends to ±00. Hence the similarity surface be- 
comes spacelike (timelike) for the positive (negative) 
mass case. We will see in Paper II that the positive- 
mass singular behavior is associated with black hole and 
white hole singularities, while the negative-mass one is 
associated with naked singularities. 

It should be noted, however, that the mass of the sin- 
gularity is not constant but is proportional to the time 
coordinate t. There is a 2-paramcter family of solutions 
belonging to this class. 



V. SUMMARY AND DISCUSSION 

In this paper, we have analyzed spherically symmetric 
self-similar spacetimes. We have shown that the metric 
of these spacetimes is conformally static, with the con- 
formal structure being determined by the form of the 
velocity function V. The homothetic Killing vector here 
plays an analogous role to the Killing vector in the static 
case. We have expressed the Einstein field equations in 
the comoving approach as ordinary differential equations 
and analytically investigated the possible asymptotic be- 
havior of perfect fluid solutions with p = (7 — l)/i and 
< 7 < 2/3. In Paper II, we focus on solutions which are 
asymptotic to the flat Friedmann universe at large dis- 
tances because we are interested in solutions embedded 
in an accelerating Friedmann background [l3| . 

These solutions can be understood in terms of the 
"complete classification" provided by Carr and Coley [3l ■ 
However, this classification mainly applies for perfect flu- 
ids with positive pressure and we have extended it to 
< 7 < 2/3. We find there are eight possible asymptotic 
behaviors at small, intermediate and large distances in 
this class of solutions. However, it should be stressed that 
there are other asymptotic behaviors, such as Kantowski- 
Sachs asymptotes at z ^ ±0, and these are discussed 
elsewhere |27| . 

Tables U and [TTl summarize the results. Table U lists 
two exact solutions: Friedmann (F) and Kantowski- 
Sachs (KS). Table [iTl lists eight asymptotic behav- 
iors: asymptotically Friedmann (F) and asymptoti- 
cally quasi- Friedmann (QF) for z s- 0, asymptoti- 
cally quasi- Friedmann (QF), asymptotically quasi-static 
(QS), asymptotically quasi-Kantowski-Sachs (QKS) and 
asymptotically constant-velocity (CV) for z 00, 
asymptotically positive-mass singular (PMS) and asymp- 
totically negative-mass singular (NMS) for z z*. 

There are three key differences from the positive pres- 
sure case. First, there is no sonic point in the negative- 
pressure case and hence no additional ambiguity associ- 
ated with the sonic point. Second, these solutions can be 
properly asymptotic to the flat Friedmann model, in the 
sense that there is no solid angle deficit. This contrasts 
with the positive pressure case, where the solutions are 
only asymptotically quasi-Friedmann. Third, the exis- 
tence of exact and asymptotically Kantowski-Sachs so- 
lutions, which can be extended beyond the timelike hy- 
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TABLE I: Self-similar exact solutions for < 7 < 2/3. The 
form of S, W and V , the number of parameters, the analytic 
continuation and the causal structure are shown. There is no 
arbitrary parameter in these solutions. 



Name 


S 


w 


V 


Continuation 


Structure 


F 


_ '2 
Z 3t 


z^ 


Z 3", 


None 


Fig.H] 


KS 


Z'^ 


z^ 


3-1 

2 1 


r — ±00 


Fig. [2] 
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persurface z = ±00 along the spacelike direction, leads 
to physically interesting solutions embedded in a Fried- 
mann background. The analytical results of this paper 
play a crucial role in Paper II, where we investigate the 
new solutions numerically and interpret them physically. 

Finally, it should be emphasized that the matter model 
considered here suffers from dynamical instability for 
small-scale perturbation. For a perfect fluid with equa- 
tion of state p = (7 — 1)^*1 the derivative dp/djj, = 7—1 
is negative for 7 < 1, which violates the condition for 
dynamically stable equilibrium. We can also discuss 
this instability in terms of linear perturbation analysis. 
Consider a density perturbation with the dependence 
^i{ut-k-x) lincaj- order in flat spacetime. Then the 
dispersion relation = (7 — l)k^ is obtained. Since the 
wave- number k is real, oj ~ — j\k\ for 7 < 1, the 

lower sign corresponding to exponential growth of the 
perturbation. Since this instability exists for arbitrar- 
ily short wavelengths, the system is unstable even in the 
presence of gravity, except for the case of a cosmologi- 
cal constant (7 = 0). This means that the matter fields 
discussed in this paper can only be valid (i.e. dynami- 
cally stable) at larger scales. In fact, there are impor- 
tant physical models which have an effectively negative 
pressure through coarse-graining but which are still ab- 
solutely stable at small scales. For example, the effective 
equations of state for domain wall networks and string 
networks correspond to 7 = 1/3 and 7 = 2/3, respec- 
tively [28j . 

It should also be noted that this small-scale dynamical 
instability does not appear for a quintessence field (i.e. 
a scalar field with a flat potential). For example, sup- 
pose the scalar potential is given by = Vb exp(V87rA(/i) 
(A > 0). This can accelerate the Friedmann universe for 
< 2. The dispersion relation for the scalar field pertur- 
bation in flat spacetime becomes uj'^ = k'^ + 8'K\'^V{(f)o), 
where 0o(t, a;) denotes the background scalar field and 
the short wavelength approximation is adopted. Thus, 
this scalar field is stable against short wavelength per- 
turbations. 
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APPENDIX A: COMPARISON WITH 
CARR-COLEY ANALYSIS. 

Carr and Coley Q have also attempted to classify 
perfect-fluid spherically symmetric self-similar solutions 
with < 7 < 1, although they have not investigated 
the solutions numerically (as we do in Paper II). Their 
analysis essentially extends the CC classification Q of 
solutions with 1 < 7 < 2 and, in principle, includes the 
< 7 < 2/3 case discussed here. However, a comparison 
with our results requires care because CC use a different 
system of notation from us. This appendix gives the re- 
lationship between the two systems and points out that 
the analysis in ref. 0] is not complete. 

CC use the same comoving coordinates (r, t, 0, 0) as 
us and the same similarity variable z = r/t, although a 
dot rather than a prime denotes d/dhiz. They also use 
the same symbols for the velocity function V , pressure p, 
density ^ and areal radial distance R = rS. Otherwise 
the mapping from our variables to theirs is as follows: 



P 



$ ^ i^, * ^ A, M ^ MS, (Al) 
^8ttGP/S^, W ^8ttGW/S^. (A2) 



They also express the field equations as ODEs in the 
quantities {S, 5", W) but they introduce the variable 



X = (47rAir2)-(i+^)/(2+7) 



(A3) 



instead of W. Note that the ODEs are expressed in terms 
of {S, M, W) in ref. 



HI 



CC correctly identify Friedmann and Kantowski-Sachs 
as the only exact power-law self-similar solutions and 
they look for more general solutions which are asymp- 
totic to these. They also point out that there is no exact 
static solutions for < 7 < 1 but they do not notice 
that there are still asymptotically quasi-static solutions 
and this is the main reason why their classification is 
incomplete. Although CC appreciate the importance of 
extending some solutions from z > to z < 0, they al- 
ways take the variable r to be positive, whereas we allow 
it to be negative. 

Their analysis of the asymptotic solutions at large and 
small distances is closely related to ours and explicitly 
includes the asymptotically Friedmann and asymptoti- 
cally Kantowski-Sachs solutions for < 7 < 2/3. How- 
ever, they do not distinguish between the asymptoti- 
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TABLE II: The possible asymptotic behaviors of asymptotically Friedmann self-similar solutions for < 7 < 2/3. The limiting 
value of z and the form of S, W and V , the number of parameters, the analytic continuation, the causal structure and the 
physical distance are shown. For the physical distance, the limiting values of zS — R/t are shown. The values of p+ = ^+(7) 
and Voo ~ 14o(7) are given in the text. 



NcLin6 
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y 


^param 


Extension 


Structure 


Distance 


F 




'2 




z^~^ 


\ 




opactjiij^e 




QF 


±0 


'2 

Z ^1 


Z^ 


z 37 


1 


None 


Spacelike 


00 


QF 


±00 


2 

Z ^1 


z" 


Z 3t 


1 


None 


Timelike 





QS 


±00 


const 


const 


2--, 
Z 1 


2 


t = ±0 


Spacelike 


00 


QKS 


±00 


z-^ 


2^ 


Z 1 


2 


r = ±co 


Timelike 


Intermediate 


CV 


±00 


z^+ 


Z 


Voo 


1 


None 


Timelike 


00 


PMS 


Zt 


2--, 

\z — Zt \ 


\z — Z* ~ 8-57 


1 


2 


None 


Spacelike 





NMS 




2--, 

\z — Z» \ ""-^T 


2-, 

\z — z»\ 8-57 


1 


2 


None 


Timelike 






cally Friedmann and asymptotically quasi- Friedmann so- 
lutions. Their description of the asymptotic-Kantowski- 
Sachs solutions is also somewhat different: because r is 
taken to be positive, they do not extend this solution 
from z > to z < 0, so they do not notice these solu- 
tions are related to wormholes (cf. Paper II). 

CC derive the asymptotically constant-velocity solu- 
tions, as well as deriving and plotting Eq. (|4.48p . but 
they describe these as "asymptotically Minkowski" and 
do fully understand their physical interpretation. They 
also identify the asymptotic solutions which represent a 
central singularity (i.e. in which 1^ ^ 1 at finite z) but 
they only include them for 2/3 < 7 < 1 and omit them 



for < 7 < 2/3. They also make no distinction between 
the positive and negative mass singularities. They claim 
to find another small-distance solution for 2/3 < 7 < 1, 
which they label "Y", but we do not find this and it is 
probably non-physical. 

Analogue of Tables I and II appears in ref. 0] but 
no analysis of the global structure of the solutions is in- 
cluded, which is another reason why the wormhole con- 
nection is missed. Carr and Gundlach [2^ have analyzed 
the global structure of some of these solutions but there 
is a mistake in the asymptotically Kantowski-Sachs solu- 
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